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The existing approximate solutions of this problem are based on 
the use of predetermined velocity profiles in transverse sections of 
the flow 1 ’ 2 . The present article proves that a solution can also be 
obtained by strict integration of equations dealing with the boundary 
layer. 


1. Integral Relations and Posing of the Boundary Value Problem 

We shall use a cylindrical system of coordinates z, R, and cp, 
and shall consider the jet emerging in the direction of the Oz axis- -the 
axis of symmetry- -from a round aperture with a radius Rj into a co- 
axial flow of the same fluid, uniform and one-dimensional at infinity, 
and possessing a velocity Uq, parallel to the Oz axis. Both in the exit 
cross section of the jet and in the cross sections of its potential core, 
we shall consider the velocity as being the same- -parallel to the Oz 
axis and equal to Uj. The region of turbulent mixing is separated 
from the jet core by surface A A' B B' C, and from the external 
potential-flow by surface A A' D D' H H' (see figure, below). The 
radii of these surfaces of revolution are designated, respectively, R 0 
and R 2 . The region of turbulent mixing, as usual, shall be divided into 
an initial sector, which extends from section A A 1 to the end of the 
core, and the main section, which lies beyond the end point C of the 
flow core. 


H 



1 H. B. Squire and I. Trouncer, ROUND JETS IN A GENERAL 
STREAM, Reports and Memoranda, A. R. C . , 1944, No. 1974. 

2 G. N. Abramovich, THEORY OF TURBULENT FLOWS (Teoriya 
turbulentnykh struy), Fizmatgiz , I960. 
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If we overlook changes in the jet pressure and consider the mean 
flow as established, we can obtain relationship 3 ’ 4 

R 2 r 2 

2 77 J* pV z 2 RdR - 2 ttU 0 f P V z RdR = pUj (u i - U 0 )irRi , (1.1) 

0 0 

which reflects the law of the conservation of linear momentum (where 
p is the mass density of the fluid and V z is the mean value of the 
longitudinal component of velocity at point z, R). However, instead of 
using the mean components of total velocity, V z and Vr, it is more 
convenient to use the components U and V, defined as 

U = V z - U 0 , V = Vr . (1.2) 

By substituting (1.2) in (1. 1), we get the corresponding equation 
for the conservation of linear momentum expressed in terms of ( 1 . 2 ): 

R 2 Rz 

2t rj P U 2 RdR + 2irU 0 J pURdR = pUi (Uj - U 0 ) wRi . (1.3) 

0 0 

We shall assume the following boundary conditions: a smooth 

transition of longitudinal velocity at the jet boundary to a velocity of 
a uniform, unidirectional flow, and the continuity of the flow along its 
axis (M-const)-- 

U U< M ' alLo = °' v Ir=. = °' u 

(1.4) 


r=r 2 


= 0 , 


au 

3 R R=R 2 


= 0 


2. Dimensionless Equation of Motion 


We select radius Ri of the exit cross section as the scale of 
length and the corresponding velocity Ui at the exit point as the scale 
of speed. Passing to dimensionless variables, we get: 


x 
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Ri 


etc; 


( 2 . 1 ) 


3 Squire and Trouncer, loc. cit. 

4 Abramovich, loc. cit. 
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( 2 . 2 ) 


The quantities Jj and Qj, having dimensions of momentum (flux) 
and mass (flux), are nondimensionalized using J 0 and Qq, the momentum 
(flux) and mass (flux) at the exit cross section. 


Thus: Jq = pU 2 jTrRi 2 , 


Qo = pUnrRi , 


(2. 3) 


and the corresponding dimensionless quantities are 


Ji 

p Ui 2 irRx 2 ’ 




Qi 

pU ittRi 2 


(2.4) 


2 Z 

Dividing both parts of equation (1.3) by pU i ttRi and using (2. 1) 
and (2. 2), we obtain the following statement of conservation of 
momentum: 


r2 r2 

2j* u 2 rdr + 2 \J urdr = 1 - \ 
0 0 



(2. 5) 


For an established mean flow of an axisymmetrical jet, ignoring 
changes in pressure and using Prandtl's formula for turbulence friction, 
we obtain the following equations 5 : 


V z 


avz 

0z 


+ 


VR 


av z 

0R 


l_ _d_ 
R 0R 


R l 


/av z \ 

\ 9R / 


av 2 , avp , v R 

3 Z 0R R 


= 0 , 


( 2 . 6 ) 


which may be restated in terms of (1.2) as 


Un + u + v 

U ° 0z 0z 0R 


JL _L 

R 0R 


H«r 


0U 0V V 

"Fz 0R R 


= 0 , 


(2.7) 


5 L. G. Loytsyanskiy, FLUID AND GAS MECHANICS (Mekhanika 
zhidkosti i gaza), Fizmatgiz , 1959. 
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where l is the mixing length in the main section of the jet, regarded 
as being proportional to R 2 , stated 6 as 

l = cR 2 (c-const), (2. 8) 


Substituting (2. 8) into (2. 7) and using (2. 1), (2. 2), and (2. 5) as well 
as taking into account that R 2 is a function of z only, we get the follow- 
ing equations for the main section of the jet: 


X 


du 

3x 


+ u 


3u 

3x 


+ 





(2.9) 


9u 9p i) 

3x 0r r 


0 . 


( 2 . 10 ) 


3. Criterion of Similarity 

The general character of the mean flow at large distances 
from the exit section of the jet ought to be determined by the quantities 
J, Uo, and p from which, however, it is impossible to construct a 
dimensionless parameter (J, the total momentum of the jet). From 
(1. 1) we have J = pUj (Uj - U 0 ) irRi 2 . From this and (2. 3) and (2. 5) 
it follows that 

J = J 0 (1 - X). (3. 1) 

Therefore, the assignment of a value to J is reduced to assigning 
values to J 0 and to X, and- -as stated previously- -the determining 
factors here are X, J 0) U 0 , and p, from which it is impossible to 
deduce any dimensionless parameter other than X itself. 

For this reason, the criterion of similarity of mean flows in two 
turbulent jets, 1 and 2, is the existence of the following condition: 



One of the conclusions which can be drawn from (3. 2) is the well- 
known fact that all turbulent jets emerging into a stationary fluid are 
dynamically similar (for all such jets, X = 0) . As is known, the laws 
governing the motion of the fluid (in the mean) in turbxilent jets 
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Squire and T r ounce r, 


loc .cit. 
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emerging into a homogenous collinear flow 7> 8 , depend primarily upon the 
value of \. 


4. Form of the Bounding Surface of the Turbulence Area 

Let us consider that the equation of boundary surface of a 
turbulence area can be given in the form of a series 


r 2 = qx^ + . . . . (4. 1) 

(On the right side of this expression is entered only the leading term of 
the expansion of r 2 in powers of x, where a is a dimensionless coeffi- 
cient. ) 


Introducing the quantities 
r 2 


1*2 


( u > = 77 / urdr > <X) = 77 / u 2 rdr | 2 <u) 2 


we may write (2. 5) in this form- 


(4.2) 



+ 


<«>- 


1 - X 

i z *z 2 


= 0 


(4. 3) 


where (u) and (u 2 ) are clearly mean quantities associated with flux 
of the fluid (mass) and flux of momentum. [In (4. 3) the coefficient 
which, according to (4. 2), takes into account the difference between 

(u) and (u 2 ) 2 , is close to 1.] From the quadratic equation (4. 3) it 
follows that 



l 



(4.4) 


Since the turbulence area spreads as the distance from the exit 
section increases, it follows that when \ > 0 and when the distance 
from the exit point is also sufficiently great, the second summand 
under the radical in (4.4) becomes much smaller than the first one. 


7 Squire and T r ounce r, loc. cit. 

8 Abramovich, loc. cit. 

9 Here and in all that follows, ( ) designates means. 
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— 2 

Expanding the radical in a series of powers of r 2 , we get: 


<»>=TT + 


r 2 


az 
r 2 4 


+ . . . , 


a2 = a 2 (\) 


(4.5) 


From (4. 5) there follows an estimate of the order of magnitude of u 
when x is large: 


( u) ~ r 2 ” 2 . (4.6) 

In estimating the mean value of v with the aid of (2. 10), we find 

( u) ~ ( u) r 2 / x . (4.7) 

In estimating orders of magnitude for the various terms, we take 10 

a 


u ~ 


/ V , v 3 1 0 111 

\ / \ / 9x x 3 r r 2 r r 2 


And in estimating the order of magnitude of various individual summands 
of the left-hand side of the equation of motion (2.9) and using (4. 7), we get 


3u (u ) 
3x X 


u 


3u (uy 


3x 


u 


3u (u)‘ 


3 r 


(4.8) 


From (4. 6) and from the fact that r 2 increases with x, it follows 
that when x is sufficiently large, <(u) is small, hence, 


(u) 2 C (u) . 

The expressions (4. 8) and (4. 9) lead to the estimate 


(4.9) 


3u 3 U ^ ^ 3 U 

u -r + l) ~Z < x , 

0x 0r 3 x 


which shows- -along with (4. 8) --that the left-hand side of the equation 
(2.9) is on the order of x “ 1 ^u^ . Both summands of the right-hand 
side of equation (2.9) are of the same order, r 2 -1 2 . Comparing 

the orders of magnitude of both sides, we get 

-l 




r 2 x 


(4. 10) 


10 


L. D. Landau and E. M. Livshits, MECHANICS OF CON- 
TINUOUS MEDIA (Mekhanika sploshnykh sred), Gostekhizdat , 1954. 
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And from (4. 10) --as well as (4, 6 ) and (4. 1 ) — it follows that 

r 2 = ax ^3 + . . . . (4. 1 1 ) 

Substituting (4. 11 ) into (4. 5), we find 

( u) = b lX “ 2 /3 + b 2 x' 4 /3 + . . . , bi = bi(\). (4.12) 

From this it can be assumed that all dimensionless functions 
defining jet flow are representable in series of powers of x“ 2 / 3 . (This 
is, in fact, later confirmed.) Hence, 

r 2 = q(\) [x 1 ^ + (3 (\)x-V3 + Y (\) x " 1 + . . .]. (4. 13) 

The initial series (4. 5) is valid when X is sufficiently greater than 
0 and when the values of x are sufficiently large. For this reason, an 
adequate correlation of the resulting series can be expected only 
within the range of the values of X and x as indicated above. When X is 
close to zero, it is difficult to estimate the order of magnitude of the 
summands under the radical in (4.4) or to determine the type of series 
for (u). As a result, one cannot successfully solve this problem when 
X is close to zero. For the case of an axisymmetrical jet emerging 
into a stationary fluid (i. e. , for the case when X = 0), there exists the 
Tolmin solution 11 . Index k in equation (4. 1 ), defining the contour of 
the boundary surface, in this case equals 1. It is obvious that in the 
range of small values of X, the power index k in equation (4. 1 ) varies 
from k = 1 when X = 0 to k = V 3 , for levels of X for which the solution 
given here already applies. 


5. Expansion in Series of the Equations of Motion and 
of Those Defining Boundary Conditions 

Equations (2. 10), (4. 12), and (4. 13) show that dimensionless 
flow function ijj, which is connected with the components of velocity u 
and v by the equations 


u 


I 

r 0r ' 


V 


1 3 ^ 

7" 0X - 


(5. 1) 


can be found in the form of an expansion of 
4 ^ ~ 4*1 4 j 2 x 2 ^ 3 + i)i3x V 3 4 • • • • 


( 5 . 2 ) 


n Abramovich, loc. cit. 
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Here the functions 4*1 depend only on the variable r\, which is 
determined by the equation 

r r 

^ r 2 a |x x /3 + (3x + yx~ 1 + . . . | 

= JL I X -V3 . (3x-i + (p2 - v) x' 5 / 3 + ...I . (5.3) 

a I 1 

With the aid of equations 

4i* = 4i > 4*2'*' = 42 c 4 ^ 3 > 4 3 * = 43 c 8 / 3 • • • , (5.4) 


a* = ac ' 2 / 3 , 0 * = pc 4 ^ , -y* = yc 6 / 3 ... 

(5.5) 

x*=xc 2 , 0 (-)/ 3 x = c 2 0 (. . ,)/ 8 x* 

(5.6) 

it is possible to convert to new dimensionless functions 4 i* as well as 
to coefficients a*, ( 3 *, y* . . . and the longitudinal coordinate x*, thus 
excluding the structural constant c from all equations. 

By substituting (5.4), (5.5), and (5 . 6 ) in (4. 1 3) , (5.2), 
we find 12 : 

and (5. 3) , 

r 2 - cl | x ^3 + px + yx^“ l + ... | , 

(5.7) 

4 = 4 i + 4 2 x " 2 /3 + 4 ^ 3 x” 4 ^ 3 + . . . , 

(5. 8 ) 

q = a ” 1 r [x“ - (3x 1 + (|3 2 - y) x” 5 /3 -f- . . . | . 

(5.9) 

We now use (5. 6 ) and (5. 9) to obtain the following for our operations of 
differentiation: 

4 - = - | x ' 1/3 - P x ~ 1 + (P 2 - y) x_5/3 + • • -I 4 - , 

3 r a 1 n 1 drj 

(5. 10 ) 


-V 3 

+ • * • ] 11 } * 

(5. 11) 


lz For the sake of simplicity, from now on we shall omit the 
asterisk from 4*1 > 42 > 43 > a > \» and x. These values will be 
presumed to be those defined by equations (5. 4), (5. 5), and (5. 6). 
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With the aid of (5. 2) and (5. 8) through (5. 11), we obtain the 
necessary series 

U = — |$ lX - 2 /3 + [® 2 - 2p$!']x- 4 /3 + [$ 3 - 2p$ 2 

a 2 

+ (3p2 - 2y)$J X" 2 + ... | , (5. 12) 

= _L {^'x- 1 + [*2« - 3p$!']x- 5 /3 + [*3«- 3P^2 ’ 

+ (6|3 2 - 3y) *!'] x“% +...}, (5. 13) 

|^ = | $i"x- 4 /3 + [® 2 " - 4p$i"] x - 2 + [® 3 » - 4(3$ 2 " 

+ (10p 2 - 4y) «!"] x" 9 /3 +...(, (5. 14) 

= -ZL - 2®i] x' 5 /3 + [-®2'ti - 4$ 2 + 5p®i'r| + 8p®i]x" 7 /3 


+ [- rj - 6^3 + Sp^z’ri + 12(3^2 “ (14p 2 - 7*y) &\ l r\ 

- ( 1 8p 2 - 12y) ® i] x - 3 +...}, (5. 15) 

U = la { ^ 1,x " 4/3 + [ 4*2' " 3 P4>i* + x ‘ 2 

+ 1 1 ^ 3 ' - pi|i 2 ' +( 5 p 2 - 5 -y) 4 »l' + — ^ j X “ 8 /3 + . . . | ( 5 . 16 ) 

wherein, for purposes of brevity, are introduced- -alongside the 
functions i(ji--also the functions 


$1 = 4*i* /t| 


(5. 17) 


The prime mark (') used in (5. 13) through (5. 17) and everywhere 
hereafter indicates the derivatives for r|. 

By substituting (5.7), (5.9), and (5. 12) through (5. 16) into the 
equation of motion (2. 9) and then adjusting the coefficients for equal 
powers of x in both the left-hand and the right-hand parts of this 


9 


equation, we obtain a group of ordinary differential equations for 
functions i}ji--as well as for functions , which are connected with 
the former by (5. 17). 

Thus, in adjusting the coefficients for , we get the following 

equation for the function 

V 3 \a 3 [c&jTi 2 ]' = [$i' 2 q] - . (5.18) 

After adjusting the coefficients for x" 7 /3 and using equation (5. 18), 
we can obtain a differential equation for function i |i 2 , 


yj- [<V$ 2 , n]' “ y [^ 2 'ri + 2 i \ j 2 ]' 


2 

To? 


4V_ 2 

*1 



(5. 19) 


while adjusting the coefficients for x“ and using (5. 18) and (5, 19) leads 
to a differential equation for function 1^3 , 


[Z^'Sj'q]' 


\a 3 

3 


[4/ 3 'q + 4i|; 3 ] 1 


\a 3 p 

3 


[ 4 , 2 - q]' + 2\a 3 P4i 2 ' 


- (^2 ,2 r,]» 


^ IMJ' 



(p 2 + ^ ) [vjji'q]' 


+ 


8a 

T~ 


4j 2 '$i 


2ap 4»i' 2 
3 q 


(5. 20) 


From (5.3) it follows that the flow axis on boundary rj has the 
following values: 


h |r=o ~ 0, h Ir = r 2 “ 1 


(5.21) 


Converting the foregoing to dimensionless values with the aid of 
(2.1) and (2.2), for conditions of (1.4), and using (5.21), we get the 
following boundary conditions: 


U I =0 <M, 


0U 

9 r 


9 , U | q=o 0. u |q= i — 0, 


0U 

0r~ 


| = 0 (M = const) 


n=l 


(5. 22) 
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We now substitute (5. 12) into (2.5) and with the use of (5.9), (5. 17), 
and (5. 21) find that the left-hand part of (2. 5) contains terms with x° , 
x~ 2 / 3, x~ 4 /3 . . . , while the right-hand side contains one term with x° . 
Adjusting the coefficients for equal powers of x on both sides of the 
equation, we get 


2\J 4V dr, = 1 - X, J 


+ w 


dr| = 0 


/[ 


- 2W 

a 2 rj 


+ 


drj - 0 . 


(5.23) 


Substituting series (5. 12), (5. 13), and (5. 16) into (5. 22), using 
(5. 17) and adjusting the coefficients for similar powers of x in the 
left-hand and right-hand sides of the resulting expressions, we find 
the boundary conditions for functions and i to be as follows: 


Irpo <M = const, ®i , | T , =0 =0, 4 ji , | T) =o =0, 4ji'| T)= i = 0, 

^i" | r|=i = 0; $2 | r|=o <M, $ 2 |r|=o = 0, Vz 1 1^=0 = 0, 


r\ 1 4^2 | n = o = 

0, V 2 * | r|=i = 0, 4J2 ,, |t 1 =i = 0; 

$3 |r|=0 <m > 

(5.24) 

O 

(1 

o 

II 

P" 

4 J 3 l |r|=o - 0» r,- 1 V3|r,=0 = °» 

4^3* lr|= 1 = 


+ 3'%=, = 0. 





6. First Approximation 


Integrating (5. 18), we get 

V 3 \a 3 T) 2 = ^i' 2 ^ + Ci . (6. 1) 

By assuming, in (6. 1) that rj = 0 and using (5. 24), we find that 
the integration constant Cj = 0, whereupon (6. 1) takes this form: 

$x' 2 = VsXa^m . (6.2) 
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In the cross sections of the flow 3u/0r ^ 0 (assuming \ <) and 
a > 0. Therefore, it follows from (5. 13) that ^ 0. Taking this 
into account, and extracting the root from both parts of (6. 2), we get: 


®i'= -V V 3 \a 3 T|$i . (6. 3) 

Integrating (6, 3) for rj, we find = C 2 -^VzyXd 3 ^ 3 . From this 

it follows, by using (5. 17) that 

4*i' = = 9 [c 2 - VV 2 7 \a 3 -r| 3 | or |l - r| V 2 ] 2 . (6.4) 

Integration constant C 2 was determined under condition 

Constant a is determined from the integration condition (5. 23) 
into which we substitute (6.4) and on computing the integral find 

a = [105(1 - X) / X 2 ] V 3 . (6. 5) 

With the aid of (6. 4) it is easy to see that all remaining conditions 
of (5. 24) are thereupon automatically satisfied. 


7. Second Approximation 

With the aid of (6.4), we get 


J 


^.dT,= " 


2- 729 


16 

7 


1 - ~ r) V 2 + — g- r) 


12 „ 3 


16 9 / , 1 6 

ry/2 + — ri 

13 1 4 n 


+ const. 


(7. 1) 


By integrating (5, 19) and utilizing (7. 1), we get the following 
equation: 


y [$2'$i ,T l] - -y [^ 2 'ri + 24 j 2 ] = ipi'ri + y 


.24 
X a 2 


87 


1 1 ^ 3 fz 

1 “ — 9 Z 2 


12 3 16 9 / 1 6 

— n - — n/2 4 — ^ 


+ d. 


(7. 2) 


If we assume in (7. 2) that r) = 0 and use (5. 24), we see that the 
integration constant D! = 0. Taking this into account and using (5. 17) 
and (6.4), we obtain from (7.2) 


12 



(i - T) 3 /2 ) (l|j 2 " T) - 4»2 ') + 3 /z T) ^2 (4j 2 't) + 24*2 ) = - Jg Xa 3 pT] 5 / 2 (l - T| ^ ) 


- d« (’ - T "* + T ” j 

- 77 t > % + T ’l 6 ) • < 7 - 3 > 

The coefficients of the linear differential equation (7. 3) and its 
right-hand side are polynomials, the subsequent terms of which contain 
r| in a power which is greater by 3 / 2 than the preceding term. The 
structure of the polynomials is such that a general solution of equation 
(7. 3) can be sought in the form of a product of q ^ and an infinite 
series in powers of 3 / 2 . It is easy to see that the infinite series in 
the present case becomes a polynomial, and that the solution looks as 
follows: 


4^2 “ Air| 2 + A 2 r| 7// 2 + A 3 t| 5 + A 4 q 13 /2 4* A 5 q 8 . 


(7.4) 


Here, A^ represents the dimensionless constant coefficients (All 
A^ are alike in that they are equal to zero, when i > 5). Solution of 
(7.4) satisfies all boundary conditions of (5. 24) when q = 0. By 
substituting (7.4) in (7. 3) and adjusting the coefficients for the highest 
power of r\^z in both parts of the equation, we find 

As = 8-729- 11 XQ ' 5) 

Adjusting the coefficients for r| 7 , r\ ll /z , r| 4 , and y\lz and utilizing, 
consecutively, (7. 5) , . . , (7. 6), we then get: 


A. = - 


A? - - 


4-31 


\q 4 , 


729-121- 13 
729-7-121 X “ 81 X ’ 


58 4 1 3 

As “ 243-5-121 Xa + 81 Xa P 


(7. 6) 


A, = 


.4 j. 1 \ _ 3 ^ 


729-121 Xa + 81 Xa P 


To determine constant (3, we shall now turn to the integral 
expression (5. 23) , which on the basis of ^2 | q = o and of (7. 1) , assumes 
this form: 


h=i 



d4 = - 


Xa 4 

120-7-13 


(7.7) 
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From (7. 4) and (7. 7) it follows: 


Ai + A 2 + . . . + Ac = - 


\a 


(7. 8) 


120-7-13 - 

Substituting expressions (7. 5) and (7. 6) into (7. 8), we get 

(3 = - 3 /izi a = -0. 024793 a . (7. 9) 

And by using (7. 5), (7. 6), and (7. 9), we find 

Ai = - 2. 9475 • 1 0” 4 \a 4 , A 2 = 2. 0730-10' 4 \a 4 , A 3 = 8. 84 2 6 • 1 0" 5 \a. 4 

(7. 10) 

A 4 = - 1. 0813- 1 0 -4 \a 4 , A 5 = 1. 5588-10 _5 \a 4 . 

(7.4) and (7. 10) will prove to us that the remaining boundary 
conditions (5. 24) are satisfied automatically. 


8. Third Approximation 


Equation (5. 20) is integrated in the same way as equation 
(5. 19) for the second approximation. The solution appears as a 
polynomial, thus 


4^3 - Bir| 2 4 B 2 r) 7 /2 4 B 3 t| 5 4 B 4 r| 13 /2 4 B 5 r| 8 4 B^r) 1 ^ 4 B 7 r) U 
wherein the coefficients have the following values: 


( 8 . 1 ) 


Bj = - 4. 8905 -10"° Xq 5 , B 2 = 5. 8683- 10" 6 \a 5 , 


B 3 = - 1. 9689 1 0“ 6 \a 5 , B 4 = - 5.9013-10 -7 \a 5 , 


-7c _5 


B s = 5. 0840* 10” 7 \a 5 


B 7 = 1. 0291- 1 0“ 8 \a 5 , 


B 6 = - 1. 3244- 10” 7 \a 5 , 


( 8 . 2 ) 


From the integral expression (5.23), we determine the constant y: 
y = - 2 . 2010 - 10 - 4 q 2 (g>3) 

We also determine another value which will be needed later: 

i|i 3 |_ = i = Bj + B 2 + . . . + B 7 = - l.l950’10' 6 \ a 5 . (8.4) 
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9. Basic Characteristics of the Flow 


The excess expenditure of fluid Q through the cross section of 
the flow is determined by the equation: 



(9. 1) 


Dividing (9. 1) by pUi u Ri 2 and using (2.1), (2.2), (2.4), and (5. 1) , 
we obtain an expression defining this dimensionless excess expenditure 
as follows: 

R 2 R-2 

q = 2 J ur dr = 2 J* dr = 2 [l(j (x, r 2 ) - 4 ( x > 0)] • (9. 2) 

o o 


If we assume that 4i| t|=o = 9, which can be done because 4i can be 
accurately determined up to the point of an arbitrary constant (see 
paragraph 6), and if we also use the conditions 42 |^=o = 43| q=o = 9> 
we can get from (9. 2), (5. 8), and (5. 21) 


q = 2 [411 1^=1 + 42 1 r|=i x “ 2 /3 + 4*3 1 7-1=1 x_4 /3 + • • • 


(9. 3) 


while it follows from (5.23) and the condition 4l|r|=o = 9 that 
1 - \ 


Hn=i = 


2\ 


(9.4) 


We now introduce the value X with the aid of this equation, 

X = x*a*“ 3 /2 = xc 2 a* _3/ 2 . (9-5) 

and by using (9. 4), (7.7), (8.4), (6.5), and (9. 5), we obtain from 
(9.3) 


q = 


1 - \ 


0. 019231 0. 00025095 + 


X 2 /s 


X 


Vs 


(9. 6) 


With the aid of (5. 7), (6.5), (7.9), (8.3), and (9. 5), we find: 
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0. 024793 


0. 00022010 


+ . . . 


(9.7) 




Vl05 (1 - X) 


X x / 3 


xVs 


X 


By assuming rj = 0 in (5. 12) and using (5. 17), (6. 4), (6. 5), (7. 4), 
(7.9), (7.10), (8.1), (8.2), (8.3), and (9. 5), we find the exce ss 
dimensionles s axial velocity of the flow as being 


u 0 


0. 037037 0. 0012470 + 0. 000045592 

x% x-/> x 2 



(9.8) 


while if we as sume T| = 1 in (5 . 1 6) and use (5 . 24) , (6.5), (7.7), (7.9), 
(8. 4), and (9. 5), we obtain an expression for the dimensionless 
transverse component of velocity at the flow boundary u*: 


V 


c 2 \ 


0. 00018315 0. 0000093208 + 


X 


X 


% 


(9.9) 


Let us determine the profile of longitudinal velocities as defined 
by two first approximations. Retaining the first two terms in (5. 12) 
and using (5. 17), (6.4), (6.5), (7.4), (7.9), (7.10), (9.5), and (9. 8), 
we get 


u 

uo 



2 


+ - 0. 0I2258t 1 3 /2 40. 027864t] 3 - 0. 018981n 9/2 +0- 003375^ _ 1Q) 

X% 


Retaining in the right-hand side of (9. 10) the one first summand, 
we obtain a velocity profile determined by the first approximation. 
Experimentally determined profiles are close to this theoretical 
profile 13 ’ 14 . The second summand of the right-hand side of (9. 10) 
turns out to be negative for the entire interval of 0 < rj < 1 . This means 
that the velocity profiles defined in the second approximation are some- 
what "sharper" than the velocity profile obtained in the first approxi- 
mation. It must also be borne in mind that the second summand 
becomes appreciable only when the value of X is very small. Generally 
speaking, at a great distance from the exit cross section of the flow, 
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13 Squire and Trouncer, loc. cit. 
14 Abramovich, loc. cit. 


it is sufficient--in all the expansions obtained- -to retain only the first 
terms, resulting in very simple expressions. Subsequent approxi- 
mations need to be taken into account only to determine the flow profile 
at the beginning of the main section of the jet flow. 


10. Starting Point for the Longitudinal Coordinate 

In the solution as found above, one thing remained undetermined 
and that was the starting point of coordinate z. A rational choice of 
such a starting point and determination of its position with respect to 
the exit section of the flow can be made only after a solution for the 
initial section of the flow has been obtained and coupled with the solu- 
tion for the main section. 


The position of the starting point of the coordinate z with respect 
to transverse section 1, passing through the end C of the core of the 
flow (see figure), can be determined as follows. 

Let us designate the coordinate of section 1 as z 0 (corresponding 
to the dimensionless coordinates x 0 , x 0 *, X 0 ). The dimensionless 
excess axial velocity in section 1 equals 1 - X: 

uq 1 = 1 - X . (10. 1) 

Substituting X = X 0 in (9. 8) and using (10. 1), we obtain the 
equation 


0.037037 0.0012470 0.000045592 1 - X 

Xo 2 /3 (x 2 /3) 2 (xof *" = x 


( 10 . 2 ) 


which serves for defining X 0 . After determining X 0 and substituting 
X = X 0 in series (9. 6) through (9. 9), it is possible to find the basic 
flow characteristics in section 1. However, we refrain from doing so 
at the moment, because the flow parameters in the transition section 
1 can be determined more accurately after the coupling of the solutions 
for the main and the initial sections of the jet flow. 


Received for publication: February 9, 1965. 
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